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Abstract 

We perform a complete calculation of charge symmetry breaking effects for the re- 
action pn — > dft at leading order in chiral perturbation theory. A new leading- 
order operator is included. From our analysis we extract 5mff, the strong contri- 
bution to the neutron-proton mass difference. The value obtained, 5m s ^ v = (1.5 ± 
0.8 (exp.) ± 0.5 (th.)) MeV, is consistent with the result based on the Cottingham 
sum rule. This agreement provides a non-trivial test of our current understanding 
of the chiral structure of QCD. 



1. At the fundamental level of the Standard Model, isospin violation is due 
to quark mass differences as well as electromagnetic effects [T|2ll3] . Amongst 
the isospin violating effects in hadronic reactions the ones that are charge 
symmetry breaking (CSB), i.e. that emerge from an interchange of up and 
down quarks, are of particular interest. Their importance is due to the fact 
that the neutral-to-charged pion mass difference, which is almost entirely 
of electromagnetic origin and usually dominates isospin violating hadronic 
observables, does not contribute here. Therefore, the sensitivity to the quark 
mass difference ma — m u is more pronounced in observables related to CSB. 

CSB effects manifest themselves in many different physical phenomena such 
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as the mass splitting of hadronic isospin multiplets (e.g. m n ^ m p [2\ and 
M D o 7^ M D + [1]), 7/-decays (for a recent two- loop calculation, see |5J and ref- 
erences therein), the different scattering lengths of nn and pp systems after 
removing electromagnetic effects in pp scattering (see, e.g. the review article 
[6]), neutron-proton elastic scattering at intermediate energies [7J, hadronic 
mixing (e.g. p° — u [8] or ir° — n [9] mixing) and the binding-energy dif- 
ference of mirror nuclei known as Nolen-Schiffer anomaly [10]. Recently, ex- 
perimental evidence for CSB was found in reactions involving the produc- 
tion of neutral pions. At IUCF non-zero values for the dd — > an cross sec- 
tion were established [llj . At TRIUMF a forward-backward asymmetry of 
the differential cross section for pn — > dn was reported which amounts to 
Afb = [17.2 ± 8(stat.) ± 5.5(sys.)] x 10~ 4 [12j. In a charge symmetric world 
the initial pn pair would consist of identical nucleons in a pure isospin one 
state and thus an interchange of beam and target would have no observable 
impact so that the cross section should be symmetric. Thus, the apparent 
forward-backward asymmetry is due to charge symmetry breaking. 

A solid theoretical background for investigating CSB effects is provided by 
chiral perturbation theory (ChPT), the low-energy effective field theory of 
QCD [T3"fl4|ll5j . Especially, since electromagnetic and quark mass (strong) ef- 
fects typically contribute with similar strength, they can only be disentangled 
within a systematic effective field theory. ChPT has been recently extended 
to pion production reactions, i.e. to processes with a large initial momentum 
p ~ V m iv M n ~ 360 MeV, with M^tum) the pion (nucleon) mass. The proper 
way to include this scale in the power counting was presented in Ref. [IE] 
and implemented in Ref. [17], see Ref. [18] for a review article. Within this 
scheme it turned out to be possible to achieve a quite good theoretical de- 
scription of s-wave pion production in pp — > dii + at next-to-leading (NLO) 
order [19]; p-wave pion production in different channels of NN — > NNtt at 
next-to-next-to- leading (N 2 LO) order was investigated in Ref. [20]. These de- 
velopments in our understanding of isospin conserving pion production mech- 
anisms provide a very good starting point for studying isospin violation effects 
in pn — > dn° and dd — > air . First efforts in this direction were already pre- 
sented in Refs. [211221128] for the pn -> dn reaction and in Refs. [24ll231f26lT27] 
for dd —>■ arc . In this work we improve the theory for the former reaction. 

The neutron-proton mass difference is due to strong and electromagnetic in- 
teractions [2], i.e. 5m n = m n — m p = 5rrf}f + 5m e ™. As a result of the chiral 
structure of the QCD Lagrangian, the strength of the rescattering operator in 
pn — > drr depicted in Fig. Q^a) is proportional to a different combination of 
5m s ^ and 5m e ^ [2"Tf2"8"] (for related work on isospin violation in pion-nucleon 
scattering see [29]). Thus, the analysis of CSB effects in pn — > dn° should 
allow to determine the values of Smff and 5m e ™ individually. This was for the 
first time stressed and exploited in Ref. [21]. Consistency of these important 
quantities as determined from pn — > dir , where they control the strength of 
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Fig. 1. Leading order diagrams for the isospin violating s-wave amplitudes of 
pn — > (Itt . Solid (dashed) lines denote nucleons (pions). Diagram (a) corresponds 
to isospin violation in the 7riV scattering vertex explicitly whereas diagram (b) in- 
dicates an isospin-violating contribution due to the neutron-proton mass difference 
in conjunction with the time-dependent Weinberg- Tomozawa operator (see text for 
details). 

the isospin violating tcN scattering amplitude, with results obtained from the 
neutron-proton mass difference itself [2] employing the Cottingham sum rule 
[30] . would provide a highly non-trivial test of our current understanding of 
QCD. It was therefore quite disturbing to find that, using the values for 5m s ^ r 
and 8rrf™ from Ref . j2] , the leading order calculation of the forward-backward 
asymmetry [21J over-predicted the experimental value by about a factor of 
3 — a consistent description would call for an agreement with data within 
the theoretical uncertainty of 15% for this kind of calculatioEpl. The evalua- 
tion of certain higher order corrections performed in Ref. [21] and in a very 
recent study [22] did not change the situation noticeably — the significant 
overestimation of the data persisted. 

In this Letter we show that there is one more rescattering operator that con- 
tributes at LO. We evaluate this new LO operator and we also recalculate 
the LO contribution considered in Ref. [21] since the numerical evaluation 
in that work turned out to be incorrect |32j. The complete LO calculation 
for pn — > dir° reveals a very good agreement with the experimental data. 
Moreover, the resulting contribution is found to be proportional to 5mff only. 
Thus, a quantitative understanding of the CSB part of pn — > dn° promises an 
alternative method of extraction of this important quantity compared to that 
used in Ref. [2]. 

2. The differential cross section of the reaction pn — > dn can be expanded 
into a series of Legendre polynomials Pi(cos6). In the near-threshold region 
only the first terms are relevant 



1 It was shown in Ref. [21] that there is no NLO contribution - thus the theoretical 
uncertainty of a leading order calculation is expected to be of the order of M^jm^. 
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where 6 is the angle between the incident proton and the pion produced and 
the Ai are functions depending on the different partial wave amplitudes. Due 
to CSB effects the differential cross section is not symmetric with respect to the 
replacement 9 «-> 7r — 6 and thus Ai is non-vanishing. The forward-backward 
asymmetry is defined as 
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where we used Eq. ([I]) in the last equality. The experiment at TRIUMF was 
done very close to threshold at Ti a b = 279.5 MeV, which is equivalent to an 
excess energy of about 2 MeV or 77 = 0.17 — traditionally, the energy for pion 
production reactions is given in terms of rj, the pion momentum in units of 
the pion mass. At this energy the total cross section a = 4ttA is dominated 
by the isospin conserving s-wave pion production amplitude. At present, this 
quantity is known theoretically only up-to-and-including terms at NLO which 
implies a theoretical uncertainty of the order of 30% for the cross section 
|19j . Therefore, to minimize the uncertainty of the current study, we use the 
experimental value for a(nn — > dn^) = 252^^ ■ 77 [fib] extracted with very 
high accuracy from the lifetime of the pionic deuterium atonf 2 !. measured 
at PSI [33]. To convert this number to the reaction of interest here we may 



use isospin symmetry which gives a(pn — > dir c 
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violating effects in this relation are to be expected of natural size and thus will 
not further be considered. In addition, we include in A also the contribution 
from the p-wave production. Here we take the results of the N 2 LO calculation 
of Ref. [20J. Thus, we get in total A = lO.Ot^ ■ 77 + (47.8 ± 5.7) ■ r] 3 [/xb]. 



At the energies we consider here, the function A\ depends on the interference 
of either an isospin conserving (IC) p-wave and an isospin violating (IV) s- 
wave amplitude or of an IV p-wave with an IC s-wave. However, only the 
former piece contributes at leading order. Thus, to the order we are working, 
one can write 
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where is the deuteron mass and k n the pion momentum. Here, M 1 and 
mI°' p are the invariant amplitudes corresponding to the isospin conserving 
p-wave pion production in the 1 5o — > 3 Sip and 1 D 2 — » 3 Sip partial waves and 
M IV ' S is the corresponding amplitude for the isospin violating s-wave produc- 
tion in the 1 P\ — > 3 S\S partial wave. Thus, in the latter amplitude the isovector 



2 Note that the Coulomb corrections were already removed in the extraction of this 
quantity from pionic atoms, see, e.g., the review 
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pion is produced from an isoscalar NN pair (ij = 0). In this work we use the 
IC p-wave amplitudes of Ref. [20]. As explained in this reference, the contri- 
bution M 1 is quite uncertain and negligibly small. We therefore negelect 
its contribution in this calculation. The IV s-wave amplitude is discussed in 
detail below. 

3. Our calculations are based on the effective chiral Lagrangian [35.15] which 
reads 



£(0) = jyt 



1 I ■ \ , 9 A -, ^ 

-T ■ (TT X 7T) H — T • a ■ V7T 



iV+... , (4) 



_4F* v ' 2F W 
for the leading irN interaction terms relevant for our study. The leading 
isospin-violating terms, generated by the quark-mass difference and hard- 
photon contributions, are 

4 0) = ^ ^s^-^f rtT-ncsN-*^ N1(t 3 *»-t.™ 3 )N+. . . (5) 

with 5m n = Smff + 8rrf™. The ellipses stand for further terms which are not 
relevant here. In the equations above F n denotes the pion decay constant in 
the chiral limit, qa is the axial-vector coupling of the nucleon and N (tt) cor- 
responds to the nucleon (pion) field. More precisely, this form of the IV strong 
and electromagnetic operators is only correct at leading order and neglecting 
terms with more than two pion fields. The more generic form involves the 
low-energy constants (LECs) c 5 and f 2 (for precise definitions, see e.g. [28J). 
Also, beyond LO other strong and electromagnetic LECs will have to be taken 
into account. 

The diagrams that contribute to the amplitude M IV,S at LO are shown in 
Fig. [TJ Diagram (a) corresponds to the rescattering process in which CSB oc- 
curs explicitly in the ttN scattering vertex due to the last two terms in Eq. (jSJ). 
In diagram (b) pion rescattering proceeds via the Weinberg- Tomozawa opera- 
tor (first term in Eq. (j3J)) which produces an additional isospin violating piece 
from the mass difference of neutron and proton due to its time dependence as 
will be discussed later in this section. 

In order to understand the interplay of diagram (a) and diagram (b) of Fig. [1] 
it is sufficient to focus on the irN rescattering vertex on nucleon 1. From the 
pion production vertex on nucleon 2 we only keep the isospin structure, for 
the rest is identical for both diagrams. The relevant part of diagram (a) then 
reads 

/(., = (t» • r*> + rfrf) + ^ (r<» ■ - rfrf) . (6) 



We work at leading order in IV. Since we study an IV transition operator, 
we may therefore treat the external nucleons as identical particles — this is 
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Fig. 2. Leading-order contributions to isospin violation due to the time-dependent 
Weinberg- Tomozawa operator in the particle basis. 

not the case for the diagram (b), where the mass difference of the external 
particles plays the essential role. The evaluation of the operator Eq. (jHJ) for 
the isospin violating transition from the isospin zero initial pn state to the 
isospin zero deuteron state yields 

{If = 0|/ W |/< = 0) = 4 (5mf - . (7) 

IT 

This piece represents the complete rescattering contribution included in Refs. [211122] . 
Let us now look more closely at diagram (b) of Fig. [TJ The relevant part of 
the amplitude for this diagram can be most easily calculated in the particle 
basis as shown in Fig. [2j One gets 

(I f = 011(b) | J* = 0) = -i(/ bl +/ b2 ), (8) 

where I\ >1 and Jb 2 are the isospin coefficients corresponding to the diagrams 
(bi) and (b 2 ) of Fig. [2] and the factor —1/2 stems from the Clebsch-Gordan 
coefficients. Note that, since the WT operator involves a time derivative, the 
corresponding Feynman rule reads 

V$ T = -^-e abc r c {q + M w ) , (9) 

7T 

with a, b and c Cartesian pion indices and the four-momentum of the inter- 
mediate pion. Due to the explicit appearance of qo in Vwt, the final expression 
for diagram (b) of Fig. [1] depends on the neutron-proton mass difference. In- 
deed, the evaluation of this vertex for the diagrams (bi) and (b2) of Fig. [2] 
yields 



!y/2 ( + 5m N ) for diagram (&i), 
r-hM \ (10) 

-V2[^fL-8m N \ for diagram (b 2 ). 
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Thus, in the isospin violating contribution to Eq. (jSJ) the terms oc M n cancel 
while those oc survive. The non- vanishing isospin matrix element for the 
diagram (b) of Fig. [TJ amounts to 

(I f = 0\I {h) \I t = 0) = ^25m N . (11) 

Adding up the contributions of diagrams (a) and (b) we find that the resulting 
contribution at LO depends on the quark mass contribution to the nucleon 
mass difference only — the electromagnetic piece vanishes completely: 

(I, = 0|J W +/ {b) |/, = 0) = 6 5mf. (12) 

7T 

In comparison with the expression used previously (cf. Eq. (j7j)) the rescattering 
operator gets enhanced by about 30%, when standard values 5m% T = 2 MeV 
and 5m e ™ = -0.76 MeV [2] are used. 

An alternative method to derive the same result is by using the field-redefined 
Lagrangian as discussed in Refs [361371 38J — see also Ref. [39] where unitary 
transformations are used. In this formulation the pion and nucleon fields are 
redefined in order to eliminate the first term in the effective Lagrangian in 
Eq. (jSJ). This allows one to work with nucleons as indistinguishable particles. 
All terms in the Lagrangian are invariant under this transformation except 
the ones involving a time derivative such as the Weinberg- Tomozawa operator 
which generates an additional isospin violating ttN — > ttN vertex oc 5m^ that 
cancels exactly the electromagnetic contribution to this vertex oc 5m e ^. 

It should be stressed that also in Ref. [21] some effects from the neutron-proton 
mass difference were included, using the formalism of Ref. [23J . However, these 
effects appear effectively in the isospin violating ttNN vertex and are explicitly 
in conflict with the chiral structure of QCD. Therefore, they are very different 
from those discussed above. 

For the sake of completeness, we present here the tree-level invariant amplitude 
M t ree S corresponding to the LO calculation 

IV[S _ 12m 2 N g A tr r dn p , {f - p) ■ p 
M tTee - i ^ dm N J ^ {pi _ p)2 + M 2, 

where p and p' denote initial and final relative momenta of the two nucle- 
ons, respectively, and p = p/p. In the calculation we use F n = 92.4 MeV 
and gA = 1-32 (utilizing the Goldberger-Treiman relation). To get the full 
amplitude M IV,S which enters the observables, M t ^ s given above needs to 
be convoluted with proper iViV wave functions in the initial and final states, 
cf. Appendix A of Ref. [20J for a detailed description. Ideally, one should use 
wave functions derived in the same framework, namely ChPT. However, up 
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to now these are only available for energies below the pion production thresh- 
old [ID]. We therefore adopt the so-called hybrid approach, first introduced 
by Weinberg [H] , i.e. we use transition operators derived within effective field 
theory and convolute them with realistic NN wave functions [12]. 

Now we are in the position to discuss the results for the forward-backward 
asymmetry within the complete LO calculation. Using the values for the pa- 
rameters specified above and utilizing the NN wave functions from Ref. |42j . 
the result can be presented in the form 

^° = (11.5±3.5)xl0- 4 ^. (14) 

As discussed above, the calculation of the coefficient has a theoretical uncer- 
tainty of 15% which is doubled to provide a more conservative estimate. This 
uncertainty is included in the expression above. We now use the experimental 
result for Aft, [12] to extract Smff which yields 

6mf = (1.5 ± 0.8 (exp.) ± 0.5 (th.)) MeV , (15) 

where we added the experimental errors in quadrature. This is the final result 
of our analysis. At the present stage, the uncertainty in the determination of 
5mff is dominated by the experimental uncertainty for Ar,. 

In this context let us point out the following: Besides the additional IV contri- 
bution discussed in detail above there are other reasons why our result deviates 
from those of Refs. [2~Tf2"2] already at leading order. The numerical evaluation 
of the diagram (a) of Fig. [1] revealed that the value we obtain is significantly 
smaller than the one found in Ref. [21]. It turned out that the result of that 
work is too large by a factor of 4 due to an error [32] . The discrepancy of our 
result to that of Ref. [22] is an accumulation of various effects. First of all 
in Ref. [22] the isospin conserving s- and p-wave amplitudes are calculated 
within ChPT up to NLO. Thus, they come with individual uncertainties of 30 
% and 15 %, respectively — the uncertainty for the s-wave appears doubled 
for this amplitude, since it enters squared in A , while the p-wave amplitudes 
mainly contribute linearly to A\ — cf. Eqs. ([2]) and In contrast to this we 
take the s-wave amplitude directly from data, with a negligible uncertainty 
and for the p-wave amplitudes the results of Ref. [20], which were calculated to 
NNLO and are additionally constrained by data. Thus, combining these uncer- 
tainties with that for the CSB amplitude in quadrature, a total uncertainty of 
50 % arises for the result of Ref. [22]- In addition, the p-wave amplitude with 
the 1 S'o initial state employed in Ref. [22], which amounts to an enhancement 
of 50% in the isospin conserving p-wave amplitude in this calculation, is in 
conflict with the data for pp — > dir + , which calls for a negligible contribution 
of this partial wave [2D]. These effects together — the larger uncertainty of 
the calculation of Ref. [22] as well as the wrong p-wave amplitude — explain 
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the discrepancy between our result and that of Ref. [22J. 

In Ref. [21j also some higher order contributions were calculated, see also 
[6]. While individually sizeable, the sum of the considered corrections was 
found to contribute very little to the asymmetry. We re-evaluated these ad- 
ditional pieces and confirmed these findings qualitatively though our results 
deviate from the ones of Refs. [2T|6] quantitatively |43j. In addition, in Ref. [22] 
some CSB p-wave amplitudes were evaluated. Through an interference with 
the isospin conserving s-wave they also contribute to the forward-backward 
asymmetry discussed in this work, however, only at NNLO. It is reassuring 
that quantitatively these contributions are in line with the power counting 
estimates given above and thus support our uncertainty estimate. 

4. In this work we calculated the CSB forward-backward asymmetry for the 
reaction pn — > dir° to leading order in the chiral expansion. We showed that the 
resulting production operator is driven by that contribution to the neutron- 
proton mass difference which is coming solely from the quark mass difference, 
Smff. Using the TRIUMF measurement of the forward-backward asymmetry 
[T2"] we extracted 

5mf = 1.5 ± 0.9 MeV , (16) 
where the theoretical and experimental uncertainties are added in quadrature. 
This number is to be compared with the value for the same quantity extracted 
from the neutron-proton mass difference — employing the Cottingham sum 
rule |30j to determine the electromagnetic contribution to the mass difference 
to 5m e ^ = -0.76 ± 0.3 MeV [2] - 

5m% r = 2.0 ± 0.3 MeV . (17) 

This value is consistent with a recent determination of the same quantity 
using lattice QCD [31J, 5mf = 2.26 ± 0.57 ± 0.42 ± 0.10 MeV was found, 
where the uncertainties emerge from statistics, from the input as well as from 
the chiral extrapolation. We emphasize that the agreement of the various 
independent extractions provides a highly non-trivial and important test for 
our understanding of the chiral symmetry and the isospin breaking pattern of 
QCD, since Eq. ffT6]) is obtained from a reaction where Smff is governed by the 
strength of ttN scattering, while Eq. (|T7j) is derived from the neutron-proton 
mass difference itself. The link between these two apparently very different 
physical quantities is provided by the symmetry pattern of QCD properly 
implemented in hadronic matrix elements through chiral perturbation theory. 

At present the uncertainty in Eq. (fl6|) is dominated by the experimental error 
bars - an improvement on this side would be very important. Still, a more 
refined calculation is also called for since only then one can be confident about 
the estimated theoretical uncertainty. Work in this direction is in progress. 
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